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The effect of Gaussian fluctuations of nanopores filled with a gas of spin-carrying molecules 
(s = 1/2) on the multiple quantum (MQ) NMR dynamics is investigated at different variances 
and correlation times of the fluctuations. We show that the fluctuations smooth out the evolution 
c ^ J ■ of MQ NMR coherence intensities which rapidly oscillate as functions of time in the absence of 

CN ■ 

fluctuations. The growth and decay of the MQ coherence clusters in the fluctuating nanopore are 
fl^ ■ also investigated. 



I. INTRODUCTION 

The multiple quantum (MQ) NMR [1] is not only a powerful tool for the investigation of nuclear spin distributions 
j^ • in solids but also an effective method to study MQ NMR coherence relaxation in large systems of highly correlated 
^yi' spins [2, 3]. With this tool, one can estimate the decoherence time [4] and a possible distance of the quantum 



X 



information transfer [2]. These parameters of quantum systems arc important for quantum information processing 

[5]. 



>, 

^O , However, dipolar coupling constants in real physical spin systems do not remain time-independent because of such 

r~^ , effects as molecular motions and imperfect experimental realizations of the two-spin/two-quantum Hamiltonian [1] 

i • together with the effect of high order corrections [6] to the average Hamiltonian theory [7] . All these effects lead to 

f^ ' uncontrolled variations of the dipolar coupling constants. A fluctuating nanopore filled with a gas of spin-carrying 

, ' (s = 1/2) molecules (atoms) is a suitable model for the study of the MQ NMR dynamics with variable dipolar 

►^ ' coupling constants [8]. This study is stimulated by experimental results concerning fluctuations in nanostructures. 

For example, fluctuating nanostructures emerge in carbon nanotubes [9] driven by mechanical or electrical excitation. 

Other examples include gas vesicles [10] and nanobubbles of insoluble gas in a liquid [11]. We expect that the results 

obtained in our paper can be useful for experimental investigations of fluctuations in nanopore materials. Because of 

the nanopore fluctuations, the dipolar coupling constant, averaged over molecular motion, fluctuates as well. Since the 

characteristic time of molecular movements is much less than the spin flip-flop times, the averaged dipolar coupling 

constant is the same for all spin pairs [12, 13] even in fluctuating nanopores. Jeener has recently shown [14] that the 

mean field approach involving the distant dipolar field [15] leads to the results which are analogous to those obtained 

in [12, 13], if equilibrium fiuctuations are taken into account in the estimation of the dipolar field. The nuclear spin 

dynamics and the NMR line shape of a nuclear spin system placed in a fluctuating nanocontainer was investigated in 
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[16, 17] on the basis of the method developed in Ref.[8]. 

In the present paper we study the MQ NMR dynamics of spin-carrying (s = 1/2) molecules (atoms) of a gas in the 
fluctuating nanopores. We investigate the relaxation (decoherence) of the MQ coherence intensities using the method 
developed in Refs. [2, 3] and applying the theoretical approach explored in Refs. [8, 18, 19]. The latter allows us to 
study the MQ NMR spin dynamics in a system of hundreds of equivalent spins. In particular, it becomes possible 
to investigate the growth of the MQ coherence clusters in a fluctuating nanopore. Assuming Gaussian character of 
nanopore fluctuations we can investigate the MQ NMR dynamics at different variances and correlation times of the 
fluctuations. 

The paper is organized as follows. The theoretical approach to MQ NMR dynamics in the fluctuating nanopores is 
developed in Sec. II. The numerical simulation of the evolution of MQ NMR coherences in the system of 201 equivalent 
spins as well as the evolution of the MQ NMR coherence clusters is presented in Sec. III. We briefly summarize our 
results and discuss further perspectives in the concluding Sec. IV. 

II. INTENSITIES OF MQ NMR COHERENCES IN FLUCTUATING NANOPORES 

We consider the MQ NMR experiment with a gas of A^ spin-carrying (s = 1/2) particles in a fluctuating nanopore 
which is placed in the strong external magnetic field Hq directed along z-axis [1] . The MQ NMR experiment consists of 
four distinct periods of time (see Fig. la): preparation (r), evolution (t), mixing (r) and detection. On the preparation 
period, the spin system is irradiated by the proper multipulse sequence of radio-frequency pulses, rotating spins by 
tt/2 around the x axis (so-called tt/2 pulses), see Fig. 16, c. As a result, the dynamics of the spin system is described 
by the effective Hamiltonian -ffe// in the rotating reference frame (Fig. la). We consider two types of pulse sequences 
on the preparation period depicted in Fig.l& and c. The first one (Fig.l&) is the standard pulse sequence resulting in 
the averaged nonsecular two-spin/two-quantum Hamiltonina Hmq [1], i.e. H^ff = Hmq- The second pulse sequence 
was introduced in [2], see Fig.lc. It yields the modified expression for H^ff which is explored below. 

Our subsequent results are based on two assumptions: (i) the characteristic time of the nanopore fluctuations is 
much longer than that of molecular motion in the nanopore and (ii) the characteristic time of fluctuations is also 
much longer than the period of the pulse sequence in Fig. 16. Then the MQ NMR dynamics on the preparation period 
is governed by the following two-spin/two-quantum Hamiltonian Hmq [18]: 

Bit) 



Hmq — -- 



.((/+)2 + (/-)2), (1) 
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where I^ = Ix + Hy and I^ = I^ — ily are the total raising and lowering operators, la (a = x, y, z) is the projection 
operator of the spin angular momentum on the axis a and D{t) is the dipolar coupling constant in the fluctuating 
nanopore, 

^(i)-7'ft^(3cos2 0(t)-l). (2) 

Here 7 is the gyromagnetic ratio, V is the volume of the nanopore, / is the form-factor, depending on the shape of the 
nanopore, and 9 is the angle determining the orientation of the nanopore with respect to the external magnetic fleld 
[12, 13]. The explicit expression for the form-factor / associated with the ellipsoidal nanopore can be found in Ref. 
[13] . Eq. (2) shows that the value of the coupling constant significantly depends on the orientation of the nanopore and 
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FIG. 1: a) The basic scheme of the MQ NMR experiment. The Hamiltonians H^ff (4) and Hmq (1) govern the spin dynamics 
on the preparation and mixing periods of the MQ NMR experiment respectively. 6) The 7r/2 pulse sequence [1]. A and 
A' = 2 A + Tp are time intervals between pulses, Tp is the tt/2 pulse duration. The period of the pulse sequence is tmq- The 
preparation period is formed by A^o cycles of the above pulse sequence, thus r = Nqtmq, Heff = Hmq, see Eq.(l). c) The 
modified pulse sequence [2] on the preparation period. The period of pulse sequence is tmq + r^z, thus r — No{tmq + Tdz)- 
The Hamiltonian Heff is given by Eq.(4) 

vanishes at the so-called "magic" angle 0,„ [20] defined as cos^m, = l/VS. We also need the Hamiltonian, responsible 
for secular dipole-dipolc interactions in the fluctuating nanopore [13], 



Hdz 



^i^ii-n 



(3) 



and the expression for the Hamiltonian Heff, governing the spin dynamics on the modified preparation period (cf. 
Ref.[2]; see Fig. Ic), 



Heff = (1 -p)HMQ+pHdz, 



(4) 



where p = Tdz/{Tdz + ta/q) < 1- 



For a fluctuating nanopore, we may write 

D{t)^{D)+5D{t). (5) 

Here (D) is the average dipolar coupling constant while SD{t) is its fluctuation. Hereafter we assume the Gaussian 
fluctuations of the coupling constant D{t) [8] and characterize the Gaussian random noise, SD(t), by the first two 
moments 

{SD{t)) = 0, {SD{h)SD{t-2)) = {i6Df)Ci\t, - fal). (6) 

Here {{SD)'^) is the fluctuation variance and C{t) denotes the correlation function. In many cases [21] 

C{t) = exp(-i/Tc) (7) 

where Tc is the correlation time. 

An important property of the Hamiltonian Hmq (1) is that it commutes with the operator of the square of the 
total spin angular momentum /^ [18]. Since the operator P commutes with I^, it is suitable to study the MQ NMR 
dynamics in a nanopore using the basis of common eigenvectors of the operators I^ and I^ [6, 18]. This choice of 
the basis allows us to split the problem into a set of simpler problems for different values of I^ . This resolves the 
problem of the Hilbert space dimension growth with the increase in the number of spins. Owing to this simplification 
it becomes possible to investigate the MQ NMR dynamics in systems consisting of several hundreds of spins and to 
calculate the profiles of the intensities of MQ NMR coherences (the intensities of MQ NMR coherences versus their 
order) [6, 18]. Eq.(l) shows that the same simplification is also valid for a fluctuating nanopore. 

In order to investigate the MQ NMR dynamics, one should first find the density matrix p[t) solving the Liouville 
evolution equation [22] 



dp 

i—- 

dt 



^Dit)[iI+r + {I~r),p{t) 



(8) 



with the initial density matrix p(0) ~ ~~Mi^ + ~, — ^z) written for the high temperature approximation [22]. Here k is 
the Boltsman constant, T is the temperature and cjq = 7-ffo- Introducing the new variable 

m = j^ 1^ Dit')dt\ (9) 

one can rewrite Eq.(8) as follows: 

^^^[HMQ,PM, (10) 

H,jQ = -^{{n' + {!-)'}, (11) 

which is the usual equation (up to some notations) for the MQ NMR dynamics. The solution of Eq.(lO) was studied 
numerically in [6, 18]. A similar consideration may be given to the evolution of the density matrix under the secular 
Hamiltonian iJ^z- Following the results of Refs.[6, 18] we represent the solution of Eq.(lO), 



as the sum of the contributions pk responsible for the MQ NMR coherences of different orders k (|fc| < N) 

pM=^Pfc(^). (13) 



k 

The intensities Jk{(p) (|fc| < N) of the MQ NMR coherences were found in Ref.[18]: 

Tr{pfc(y))p_fc(y))} 
Jkiv) - ^^j2) ■ (14) 

We also consider the decay of the MQ coherence intensities in the MQ NMR experiment with the modified preparation 
period [2] (see Fig.lc) concatenating the short evolution periods Tdz under the averaged over fluctuations secular 
dipole-dipolc Hamiltonian (3) 

Ha.^^-^CMl-1') (15) 

(which is considered as a perturbation) with the evolution period tj^iq under the nonsecular average two-spin/two- 
quantum Hamiltonian Hmq (11). Then the the Hamiltonian H^ff (4) averaged over fluctuations is [2] 

H,ff{p) = {l-p)HMQ+pHd,. (16) 

In this case, the intensities of the MQ coherences Jk{(p,p) are the following [19]: 



Tr{pk{ ip,p)p-k {(p,p)} 
Tr(/2) 



MV,P)^ 7KJJ2\ ■ (17) 



Here we use the following representation for the density matrix p: 

k 

where pk is the contribution to p from the MQ coherence of the fcth order. We will use both intensities (14) and (17) 
in our calculations. 
It is suitable to write 

Jk{'p)^Jk{t + X{t)), (19) 

where 

1 /■* 



X{t) = —jJD{f)dt'. (20) 

By the definition of the Dirac delta function one has the following identity: 

/oo 
Jk{t')5{t' -t-X{t))dt'. (21) 

OO 

Using the integral representation of the ^-function, 

5{x) = ^p duje-^-\ (22) 

we transform Eq.(21) to the following one: 

1 f°° 
Jk{t + X) = — / dw Jfc(t')e~*"^* -*)e*'"^(*)(it'. (23) 

27r ./_^ 



This equation can be easily averaged over fluctuations. First of all, averaging both sides of Eq.(23) one can write: 

{Jkit + Xm = — / dwJfc(i')e-'"(* -*)(e*'^^(*))dt'. (24) 

It is known that [8, 21] 

l^^^u:X{t)^ = exp ( - uj^{{SDf)T^{t)^ (25) 

for the Gaussian fluctuations, where 

T^{t)= f {t-t')C{t')dt'. (26) 







Substituting Eqs.(25) and (26) into Eq. (24) one can obtain after simple calculations: 

(j,(t + x(t))) = „ , ,..i.,^..,. r di'Mt')oM'TT}rT^vSn7^)- (27) 



Eq.(27) will be used in Sec. Ill for study of the dependencies of the MQ NMR coherence intensities on different 
variances and correlation times of the fluctuations in the system consisting of 201 spins. Here the intensities J^ are 
defined by one of the formulas (14) or (17) and the correlation function, C(<), is taken from Eq.(7) in all calculations. 
Substituting Eq.(7) into Eq.(26) one obtains 

T\t)=tl[e-^.-,^--\\^. (28) 

Eq.(27) shows that the fluctuation effect is characterized by the function {{5Dy')T'^{t) which takes the following form 
in view of Eq.(28): 

{{5Df)T^{t) = {{6Df)Tl (e-^^ + f - l) ■ (29) 



If I/tc ^ 1 , then one has 



{{5DY)T\t) « {{5DYWt, (30) 



so that the ffuctuations are characterized by a single dimensionless parameter {{SDY)tc/{D) = ".^^A ' T(.{D) which 
is a product of the relative variance ijjyi and the dimensionless correlation time Tc{D). Another limit is I/tc ^ 1. 
Then one has 

{i6Dr)T\t) « \{iSDf)t', (31) 

so that the fluctuations are completely characterized by the relative variance {{5D)'^)/{D)'^. 



III. NUMERICAL SIMULATIONS 

The MQ NMR dynamics of spin bearing molecules in a fluctuating nanopore can be characterized by the averaged 
(over the fast molecular motion in the nanopore) dipolar coupling constant, (-D), by the variance, {{SDY) and by the 
correlation time, Tc, of the fluctuations. We use such dimensionless parameters as the relative variance A, the time i 
and the correlation time fc which are introduced as follows: 

A = ^^^, t = t{D), f, = r,{D). (32) 



A =1, Tc=0.1 
J=l, Tc=0.01 
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FIG. 2: The evolution of the zero order coherence intensity Jo(i) in the presence of the fluctuations with different relative 
variances A and correlation times fc 



First, wc consider the influence of the fluctuations on the profiles of the MQ NMR coherence intensities in the standard 
MQ NMR experiment [1] with the spin-1/2 particle system of 201 equivalent spins. One has to emphasize that the 
MQ NMR coherence intensities are quickly oscillating functions in the system of equivalent spins [6] . For this reason, 
in order to construct the profiles of the MQ NMR coherence intensities we consider the intensities Jk{t) (see Eq.(14)) 
averaged over the sufficiently long time interval 31 < t < 45.5 and find the profiles of these averaged intensities, see 
Ref.[6]. Now, taking into account the fluctuations, we observe that the fluctuations smooth out the oscillations of the 
coherence intensities. This is demonstrated in Fig. 2, where the evolution of zero-order coherence Jo(i) is represented 
for different relative variances A and correlation times fc. We see that, while reducing oscillations, the fluctuations do 
not change the mean value of the zero order coherence Jg. The same conclusion is valid for the higher order coherence 
intensities as well. Thus, one can expect that fluctuations do not change the proflles of the MQ NMR coherence 
intensities averaged over the long time interval. The latter conclusion is justifled by the direct calculations of the 
MQ NMR profiles for the system of 201 spin-1/2 particles for the fluctuations with different relative variances and 
correlation times. Since proflles are not affected by the fluctuations in the standard MQ NMR experiment, we do not 
discuss this case in more details. 

Now we refer to the experiments with a modifled preparation period, where the decay of MQ NMR coherences is 
inherent in the preparation period [2, 19]. Intensities of MQ NMR coherences are given by Eq.(17) in this case. We 
study the fluctuation effect on the evolution of the coherence clusters. By the coherence cluster we mean the set of 
coherences having intensities Jk > 0.005 [19]. It is obtained that the fluctuations effect both the maximal size of 
the coherence cluster iV™"^ = Nci{ici) and the formation time of this cluster Id- This is demonstrated in Figs. 3, 
where the evolution of the coherence cluster size is shown for two values of the parameter p {p = 0.004 in Figs.3(a, h) 
and p = 0.009 in Figs.3(c, d)) and different relative variances and correlation times of fluctuations. These flgures 
demonstrate that, in general, the evolution of the cluster size Nd depends on both the relative variance A and the 
correlation time fc. 

However, the correlation time effect disappears, if t/fc <C 1, which follows from Eq.(31). This fact is suitably 



p=0.004 



(a) 



20-, 



15- 






10- 



5- 



0- 




I ' ' ' I ' ' ' I 
10 12 14 16 



Dimensionless time, t 
p=0.009 



(4 




1 ' ' ' ' I ' ' ' ' I ' ' ' ' I 

2 3 4 _5 
Dimensionless time, t 



20 -[ 



15- 






10- 



5- 



0- 



15- 

^ 10- 



J- 



0- 



p=Q.004, Tc =1 



m 




loT^ 



4 6 S 10 12 14 16 
Dimensionless time, t 



p=0.009, Tc=l 



m 



.-J.-: 




-r-n 



2 3 4 5 

Dimensionless time, T 



FIG. 3: The evolution of the cluster size Nd under fluctuations with different relative variances A and correlation times 
fc- Cluster evolution in the spin system without fluctuations (A = 0) is shown by solid lines; (a) p — 0.004, A = f, 
fc = O.OOI, f, fOO; (b) p = 0.004, A = O.I, I, 10, f^ = I; (c) p = 0.009, A = I, r^ = O.OOI, I, 100; (d) p = 0.009, 
A = O.I, I, 10, fc = I 



reflected in Table I, where the maximal cluster sizes iV™°^ = Nci{tci) and the time moments of their formation td 
for p ~ 0.004 and p ~ 0.009, corresponding to the relative variance A = 1 and different correlation times fc are 
represented. In fact, parameters iV™°^ and td corresponding to fc = I and fc ~ 100, are essentially the same, which 
is valid for both p = 0.004 and p = 0.009. This conclusion agrees with Eq.(3I), where only the relative variance is 
represented as a parameter responsible for the fluctuations. 

On the other hand, the relative variance effects considerably on the cluster evolution no matter how big this variance 
is. In fact, the relative variance is present in all three formulas (29-31) so that it is a significant parameter of spin 
dynamics. This is demonstrated in Figs. 3(6, d) as well as in Table II, where the parameters N^'^^ and id for p — 0.004 
and p = 0.009, corresponding to the correlation time fc = I, and different relative variances A are collected. 
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TABLE I: The maximal cluster size N^"-^ and the time moment of its formation td for the relative variance A = 1 and different 
correlation times fc- The appropriate parameters of the spin system without fluctuations are the following: j^j"^"-^ — 20, 
td = 1.36 for p = 0.04 and iV^""'' = 19, id = 1.56 for p = 0.09 
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TABLE II: The maximal cluster size jV™"^ and the time moment of its formation td for the correlation time Tc = 1 and 
different relative variances A. 

IV. CONCLUSIONS 



We consider fluctuations of the dipole-dipole interaction constant D in a system of equivalent spin-1/2 particles. 
It is found that the fluctuations smooth out the evolution of the MQ NMR coherence intensities. This is the only 
fluctuation effect unless the decay process is involved. Decay of the MQ NMR coherences has been considered in the 
MQ NMR experiment with the modified preparation period [2, 19], where short periods of MQ coherence formation 
alternate with short periods of their decay. We find that the fiuctuations affect both the maximal size of the coherence 
cluster and the period of its formation. Namely, the maximal cluster size and the period of its formation decrease 
with the increase in both the relative variance A and the correlation time fc. However, the effect of the correlation 
time is reduced if ici/Tc <C 1. 

Although a nanopore is a rather specific quantum object, the numerical results obtained for it may be useful for 
study of the fluctuation effects in other objects. However, numerical simulation of the fiuctuation effects in quantum 
systems, different from nanopores, is much more complicated because of the exponential growth of the Hilbert space 
dimension with the increase in the number of the spins, involved in the quantum process. A detailed study of different 
aspects of the MQ NMR dynamics [2, 3, 18] is aimed at the increase of the power of the MQ NMR spectroscopy in 
the investigation of the properties and structures of different physical objects. Using our approach, the information 
about volumes, shapes and orientations of nanopores in materials with identical nanopores can be readily extracted. 
It is also possible to obtain information about the distribution of nanopore parameters in a system with non-identical 
nanopores [23]. In particular, information about orientations of nanopores can be deduced from the free induction 
decay (FID) using different orientations of the sample with respect to the magnetic field. In fact, the long time 
behavior of the FID is determined mainly by the nanocavities, where the dipolar coupling constants between spin- 
earring atoms (molecules) are vanishing, which corresponds to the so-called "magic" angle orientation of the nanopore. 
[20]. The obtained results can be also averaged over the orientations of nanopores. 

All numerical simulations have been performed using the resources of the Joint Supercomputer Center (JSCC) of 
Russian Academy of Sciences. The authors thank Dr. S.I.Doronin for his assistance in obtaining the numerical data. 
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